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' Kinematic fitting is a well-established tool to improve jet energy and invariant mass 

resolutions by fitting the measured values under constraints (e.g. energy conservation). 
However, in the presence of substantial ISR and Beamstrahlung, naive energy and 
(longitudinal) momentum constraints fail due to the a priori unknown amount of un- 
detected momentum carried away by coUinear photons. It is possible to take care of 
those two effects and thus obtain significantly higher mass resolutions. 
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1 Introduction 
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I ' Kinematic fitting is used to improve measm'ed quantities and/or determine unmeasured 

quantities by using known constraints of the events, e.g. momentum conservation. Already 
(— I ' at LEP this led to a significant improvement of the W mass resolution [2] . The special new 

problems at the ILC are the substantially higher amounts of initial state radiation (ISR) 
and Beamstrahlung which are due to the increase in beam energy and luminosity. 
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\0 [ 2 Benchmark Sample 

To show the problem and to test possible reme- 
dies, a sample e+e^ uddu with an integrated 
luminosity of 15 fb~^ is used. The quarks stem 
O ; mostly from W+W' and about 1% from ZZ de- 

cays. A sample of light quarks has been chosen 



Cut 


events 


% 


total sample 


52490 


100.0% 


|cos(6'jct)| < 0.989 


39 940 


76.1% 


A^track/jet > 1 


38 925 


74.2% 


Sjct > 4.5 GeV 


38 912 


74.1% 



to ensure absence of B-mesons and subsequent ^^^j^ ^. ^^^^ ^^^y^^^ j^^^^ 

leptonic decays. The sample is a full Monte- 
Carlo simulation of the LDCPrime_02Sc detector model with full reconstruction using 
$H ' Pandora particle flow [3]. The reconstructed particles are forced into four jets using the 

Durham Jetfinder [5]. In order to select well-reconstructed events, the cuts listed in Ta- 
ble [T] are apphed to these jets. Based on MC information, two subsamples are then selected 
to separately investigate events with and without a significant amount of missing energy 
-Emiss = \fs — X] £'quark, whcrc ^/s IS the nominal center of mass energy (500 GeV): 

• 'no ii^miss < 5 GeV; 15 726 events =40% of sample after first cuts 

• Eraiss > 30 GeV; cos 6*^ > 0.999, iiE^>b GeV; 9 186 events = 24% 

The remaining events (little missing energy or photons within detector acceptance) will be 
investigated later. Five constraints are applied to both subsamples: 

^1-4: four-momentum conservation: X^Pjots ~ (\/s, 0,0, 0) 

#5: equal invariant di-jet masses for those jet pairs from Wi/Zi and W2/Z2 
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Since the pairing of the jets is a priori unknown, all three possible combinations are fitted. 
The "best jet combination" is the one with the highest fit probability > and a di-jet mass 
of 50 GeV < mjj < IfO GeV. 



3 4-Jet Hypothesis used as Benchmark 




Energy of unfitted jets [GeV] Fit probability (best jet combination) 

Figure 1: Both subsamples: (a) sum of jet energies, (b) fit probability. 

Fig. [ija) shows the total visible energy for both subsamples. In the case 'no the 
constraint J2^3 — ^00 GeV can be fulfilled within some energy spread (80%), while a fit 
enforcing such a constraint has very little chance to converge in the other case (24%). This 
can be seen in Fig. [Tfb), where the fit probability for the best jet combination is shown 
(converged fits only). 

4 Modelhng ISR and Beamstrahlung 

In this first approach, possible ISR and Beamstrahlung photons are modelled as one object 
with measured parameters and errors reflecting the missing momentum spectrum. Contrary 
to other approaches (unmeasured parameters, soft constraints), no hard constraints are lost. 

The photons are parametrized by Px,Py,Pz- Since both ISR and Beamstrahlung are 
mostly emitted parallel to the beam, p^, and Py are fixed to zero and only is varied by 
the fit. The "measured" value of p^ is set to zero and the error is assumed to be Gaussian 
with a = 100 GeV (first try, since easy to implement). 

The MC simulation generates one photon each for both e"*" and e~ . However, in most 
cases only one of these two photons carries away a substantial amount of energy. Thus it is 
considered sufficient to add only one photon to the fit hypothesis as shown in the following. 

5 4-Jet-Plus-Photon vs. 4-Jet Hypothesis 
5.1 'p' Subsample 

Fig. [2] shows a comparison of the two different hypotheses on the subsample. Fig. [2ja) 
shows again the total visible energy, for 4-jet+7 the fitted photon energy is added. Now 
it is also possible for the fit to converge on the subsample (79%). The fit probability 
distribution is shown in Fig.[2ljb). 
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Figure 2: Both fit hypotheses, tested on the subsample: (a) sum of jet energies, (b) fit 
probabiUty, (c) invariant di-jet masses before and (d) after the fit. 



01200 



2 150 





1 Entries 7233 |" 




1 




























■ • - • wFi"'^ - 












1 









Missing IMC energy [GeV] 

Figure 3: Missing energy in subsample: 
fitted vs. MC. 



Fig- Etc) and (d) show the invariant 2- 
jet masses before and after the fit. In the 
4-jet hypothesis the peak of the distribution 
is shifted towards higher masses, since tlie fit 
tries to compensate the missing energy to sat- 
isfy energy conservation. Thus, too large in- 
variant masses are fitted, which makes the dis- 
tinction between W and Z more difficult. The 
4-jet-|-7 fit includes the effect from ISR and 
Beamstrahlung and diminishes the bias and 
the mass peak width significantly. From the 
correlation between fitted photon energy and 
the missing MC energy i^miss in Fig. [3] can be 



seen that the fit indeed recovers the energy of the escaping photons. 



5.2 'No Subsample 

Since, ideally, one fit should be sufficient for all events, the effects of the 4-jet+7 hypothesis 
on the 'no ^' subsample are investigated as well. Fig. [4] contains the distributions corre- 
sponding to Fig. [2] for this subsample. Fig. ID^a) shows that small amounts of energy are 
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Figure 4: Both fit hypotheses, tested on the 'no ^' subsample: (a) sum of jet energies, (b) 
fit probabihty, (c) invariant di-jet masses before and (d) after the fit. 



recovered by the 4-jet+7 fit ('no ^' contains missing energy up to 5 GeV), but the fit con- 
verges for fewer events as shown in Fig. I^b) (to be investigated). Fig.[lKc,d) ilhistrate that 
4-jet-|-7 yields slightly less resolution improvement, but eliminates the bias of the 4-jet fit. 

6 Summary and Outlook 

As the first results show, it is possible to reconstruct ISR photons and improve the di-jet 
mass resolution by about 1 GeV already with a very basic photon parametrization. Further 
improvement, especially concerning the convergence, can be achieved by a more accurate 
modelling of ISR and Beamstrahlung, better error estimation and jet energy scaling. The 
concept will be tested with other quark flavors and compared to a soft constraint approach. 
A special treatment for ISR photons measured in the detector has to be implemented. 
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